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. . . Abstract 

We have calculated evolution of a non-thermal electron population from super-thermal but weakly 
relativistic to highly relativistic energy range in clusters of galaxies. We investigate evolution of hard X-ray 
' radiation due to both bremsstrahlung and inverse Compton scattering of the cosmic microwave background 

$H , photons. The bremsstrahlung component is more significant than the inverse Compton scattering one when 

Oh 1 the momentum spectra of electron sources are steeper than ~ -P e _3 '° and vice versa in the case of Coma, 

where P e is an electron momentum. The resultant hard X-ray spectra are flatter when the bremsstrahlung 
component is dominant. When the spectral indices of the source term are in the intermediate range 
(—2.5 ~ —3.5), too much extreme ultraviolet emission is produced. Inverse Compton dominant models 
can reproduce Coma cluster results with reasonable injection rates, which are possible in cluster mergers 
^ ' and/or ambient gas accretion. 

Key words: galaxies: clusters: general — galaxies: clusters: individual (Coma) — intergalactic medium 
— radiation mechanisms: non-thermal 

O 

o 

O 



X 



INTRODUCTION 



Clusters of galaxies (CG) contain not only the thermal intracluster medium (ICM) but also non-thermal high energy 
particles in intracluster space. One of their direct evidences is existence of non-thermal synchrotron radio halos and 
Q ■ relics, which have been observed since 1970s (Wilson 1970; Giovannini, Tordi, & Feretti 1999). This indicates that 
there are non-thermal electrons with energy of ~ GeV in intracluster space. Although their origin is still not clear, 
c/3 • they are probably related with some active phenomena in CG such as cluster mergers, AGNs, AGNs' jets, and/or 
star burst in the member galaxies. Therefor, it is crucial to study non-thermal high energy particles in order to 
investigate active phenomena in CG. Furthermore, such non-thermal electrons and protons possibly have an effect on 
ICM dynamics (Rephaeli & Silk 1995; Inoue & Sasaki 2001) and Sunyaev-Zel'dovich effect (Blasi, Olinto, & Stebbins 
2000; Ensslin & Kaiser 2000). 

On the other hand, non-thermal hard X-ray radiation due to inverse Compton scattering (ICS) of cosmic microwave 
background (CMB) photons by the same electron population is expected from such CG with radio halos and relics 
(Rephaeli 1979). There have been many attempts to detect such high energy emission (e.g., Rephaeli, Ulmer, & 
Gruber 1994; Henriksen 1999; Valinia ct al. 1999). However, because of the dominant thermal X-ray emission from 
ICM, it was very difficult to detect non-thermal hard X-ray radiation. 

Recently, excess of hard X-ray radiation over thermal emission is detected from Coma (Fusco-Femiano et al. 1999; 
Rephaeli, Gruber, & Blanco 1999), A2256 (Fusco-Femiano et al. 2000), and HCG62 (Fukazawa et al. 2001) though 
their origin is still unclear. If we assume all hard X-ray radiation is due to ICS of CMB, we can estimate intracluster 
magnetic field strength or its lower limit by comparing hard X-ray flux or its upper limit with synchrotron radio 
flux, respectively (e.g., Fusco-Femiano et al. 1999; Henriksen 1999; Valinia et al. 1999). However, the magnetic 
field strength estimated through this method tends to be weaker than that determined through Faraday rotation 
measurement (Kim, Tribble & Kronbcrg 1991; Clarke, Kronberg,& Bohringer 2001). 

One possible solution of this discrepancy is that all or a significant part of hard X-ray emission is due to non-thermal 
bremsstrahlung (NTB) from weakly relativistic non-thermal electrons with energy of a few 10 or 100 keV (e.g., Kaastra 
ct al. 1999; Ensslin, Lieu, & Biermann 1999; Kempner & Sarazin 2000; Dogiel 2000; Blasi 2000). However, NTB needs 
much more non-thermal electron energy than ICS to produce a given amount of hard X-ray. In addition, dominant 
cooling process of electrons relevant to NTB hard X-ray is not bremsstrahlung losses but Coulomb interaction with 
thermal ICM, which results in energy transfer from the non-thermal electrons to the thermal ICM. Petrosian (2001) 
showed that unusual rapid heating will occur if Coma hard X-ray is due to NTB. However, calculation of NTB spectra 
is still interesting because it will possibly work well to explain future observations of other clusters. At present, hard 
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X-ray excess have been detected from only a few clusters of galaxies. Thus, it is not clear whether Coma cluster is a 
typical case or not. 

Whether ICS or NTB is dominant in hard X-ray depends on a shape of electron momentum spectra. This depends 
on acceleration and cooling processes of electrons. For instance, if turbulent acceleration in ICM is an only effective 
acceleration process, it is difficult to accelerate thermal electrons up to ~GeV because the acceleration time is too 
long for typical intracluster conditions. Therefore, almost hard X-ray is due to NTB (Blasi 2000). On the other hand, 
if we consider shock acceleration in ICM (Takizawa & Naito 2000), it is expected that accelerated electron momentum 
spectra are power-law up to the critical energy where an acceleration timescale is comparable to an ICS cooling 
timescale. For typical rich clusters, this critical energy is a few TeV (Loeb & Waxman 2000; Totani & Kitayama 
2000). Thus, both ICS and NTB can contribute to hard X-ray radiation. When electron momentum spectra are single 
power-law, Sarazin and Kcmpner (2000) shows that NTB is dominant in hard X-ray if the momentum spectral index 
is more than ~ 2.7. 

Although accelerated electron momentum spectra in shocks are single power-law, total spectra in CG do not have 
to have the same form because various cooling processes can change the spectral shape. In a higher energy range, ICS 
and synchrotron loss modify the spectral shape from that of the source. In lower energy range, on the other hand, 
Coulomb loss does. Therefor, we have to to follow time evolution of electron momentum spectra to investigate hard 
X-ray emission from CG. 

In this paper, we calculate time evolution of electron momentum spectra for spherically symmetric models of CG 
and investigate evolution of both NTB and ICS contributions in hard X-ray. It is certain that the synchrotron radio 
spectrum of Coma is known to a much higher precision than its hard X-ray spectrum. However, synchrotron cmissivity 
is proportional to the non-thermal electron density and the magnetic field energy density. Unfortunately, we do not 
have enough information about their spatial distribution. On the other hand, ICS and NTB emissivity are proportional 
to the non-thermal electron density and, the CMB energy density and the thermal ICM density, respectively. We have 
enough spatial information about both of them from radio and X-ray observations. Thus, we think calculating the 
NTB and ICS hard X-ray spectra is as useful as calculating synchrotron radio ones in order to consider physical 
property of non-thermal electrons in clusters of galaxies. The rest of this paper is organized as follows. In §[2] we 
describe our models. In ^ we present the results. In §0 we compare our models to Coma cluster and discuss their 
implications. Finally, our conclusions are presented in §pT 

2. MODELS 



2.1. Evolution of Electron Momentum Spectra 

To follow evolution of momentum spectra of a non-thermal electron population, we should solve a diffusion loss 
equation (Sarazin 1999). According to Bohm diffusion approximation, however, the spatial diffusion term is negligible 
in typical intracluster conditions within a few Gyrs (Takizawa & Naito 2000). We consider an isotropic momentum 
distribution function. It is not clear whether the electron pitch angle distribution is randomized or not in typi- 
cal intracluster conditions. However, this isotropic assumption is clearly valid when we consider Inverse Compton, 
bremsstrahlung, and Coulomb losses because Faraday rotation measurements toward clusters of galaxies suggest tan- 
gled magnetic field in intracluster space. Chaotic field lines also give rise to isotropic emission independent of the 
electron pitch angle distribution. Only synchrotron loss rate depends on the pitch angle distribution but it does not 
dominate other loss processes within our parameter range. We consider spherically symmetric ICM for simplicity. We 
will use the normalized electron momentum defined by p e = P e /(m e c). Let N(p e ,r,t)dp e be the number density of 
non-thermal electrons with momenta from p e to p e + dp e at the radius r. Thus, the equation which describes time 
evolution of N (p e ,r,t) is, 

dN{p^t) = A. [b{Pejr)N(Pe ^ t)] + Q{Pe ^ t)i (1) 

where, b(p e ,r) is the rate of normalized momentum loss for a single electron with a normalized momentum of p e at 
the radius r, and Q(p e ,r,t) gives the rate of production of new non-thermal electrons. 

We consider ICS of the CMB photons, synchrotron loss, and Coulomb loss to determine b(p e ,r). Bremsstrahlung 
loss is neglected because it does not contribute significantly in our parameter range. Thus, 

^e,r)=-f%l -(^) . (2) 



v dt J ic V dt ) sync \ dt / Coul 
Each component is numerically (sec Longair 1997), 

=4.32xlO- 4 ^-^Gyr- 1 , (3) 
\dtJ IC pl + 1 
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Fig. 1. (a)An electron momentum loss rate b(p e ) as a function of p e . (b)Electron cooling time defined by p e /b(p e ). In both cases we 
assume n e = 10 _3 cm -1 and B = IfiG. ICS loss is dominant for a higher energy region (p e S> 100), where ft(p e ) ocpj; and Pe/b(p e ) ocpj" . 
The loss rate is well approximated by Coulomb loss with an extremely relativistic limit in a middle region (p e ~ 10), where b(p e ) is 
almost constant and p e /b(p e ) oc p e . Transrelativistic effect appears in a lower energy part (p e ~ 1), where b(p e ) oc (p^ + l)/p^ and 
Pe/b(p e ) oc pl/(pl + 1) approximately. 

dp ^ = 4.10xlO- 5 ^ T ^f^V G yr- 1 , (4) 



dtj avnc • Pl + l^K^G 



Vp!+t 



±) = 3 .79xlo(,^)^(l.O + ^ + ^S>W- (5) 
dtJ Coul Vl0- 3 cm-V v\ \ 19 75 J 

where B is the intracluster magnetic field strength, and n e is the electron density of thermal ICM. We consider CMB 
photons at z=0 as soft photon density calculating inverse Compton loss rate. We use 2.7 K as CMB temperature. 

The loss rate b(p e ) and cooling time defined by p e /b(p e ) are shown in Figure [l] as a function of p e when n e = 10 _3 cm _1 
and B = l^G. ICS loss is dominant for a higher energy region {p e ^> 100), where b(p e ) ocpl and p e /b(p e ) ocp^ 1 - The 
loss rate is well approximated by Coulomb loss with an extremely relativistic limit in a middle region (p e ~ 10), where 
b(p e ) is almost constant and p e /b(p e ) oc p e . Transrelativistic effect appears in a lower energy part (p e ~ 1), where 
b(p e ) oc {pi + l)/p1 and p e /b(p e ) oc p 3 e /{p 2 e + 1) approximately. 

We assume that Q(p e ,t) ocp^^, which is expected from the standard theory of shock acceleration. We neglect the 
temporal change of the index fx, 

2.2. Emission from a non-thermal electron population 

Let l e (r) be the volume emissivity of non-thermal hard X-ray at energies from e to e + de at the radius r. We consider 
NTB and ICS of CMB photons from a non-thermal electron population, 

Ur)=lf™(r) + ll cs (r). (6) 

For NTB components, the emission is given by the integral (see Sarazin & Kempner 2000), 

ir*(r)=e I N{ Pl ,r)dpMPi)Y. n ^) da{Pl d e > Z \ (7) 
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where, v{p\) = cp\/{p\ + l) 1 / 2 is the velocity of an electron with normalized electron momentum pi, Z is the 
charge of the various thermal particles in ICM, nz{r) is the number densities of these thermal particles at r, and 
[<i<j(pi, e, Z)/de]de is the brcmsstrahlung cross section for emitting photons with energies from e to e + de for the 
collision between an electron with normalized momentum p\ and an ion with electric charge Ze. 

To calculate the cross sections, we use the Born-approximation formulas of completely unscreened nuclei (Koch & 
Motz 1959). We use the extreme-relativistic limit of Bethe- Heftier formula for the extremely relativistic (p e > 4.81) 
electrons. On the other hand, we use the transrelativistic expressions of equations (2) and (4) in Haug (1997) for 
weakly relativistic electrons (p e < 0.965) and moderately relativistic (0.965 < p e < 4.81) electrons, respectively. The 
distortion of the electron wave function by the nuclear Coulomb field is considered by multiplying the Elwert factor. 

For ICS components, the emission is (see Blumenthal & Gould 1970), 

TO—fe^/VTlVV^ (8) 



where o~t is the Thomson cross section, h is the Plank constant and 

F(x) = l + x + 2xlnx-2x 2 . (9) 
J{u) is the mean intensity at frequency v of the seed photon field. For the CMB, this is the black body function 
_ 2hv z 1 

M_ c 2 exp(hv/kT CMB )-l { ' 

The luminosity at energies from e to e + de from a whole cluster L e is obtained by spatial integration of the above- 
mentioned volume emissivity, 

L e = Lr B +Lf s , (11) 

where, 

Lf TB = J lf TB ATrr 2 dr, (12) 

L ics = f £CS 47rr 2 dr . (13) 



2.3. Model Clusters 

For the thermal ICM, we assume the isothermal beta model (Cavalicre&Fusco-Femiano 1978). The electron number 
density of ICM at the radius r is, 



n e (r) = n e ,o 



1- 



(14) 



and the temperature T is constant. In order to compare our results with Coma cluster in §Q, the parameters of 
Coma are used. We set n = 3.12 x lO^cm -3 , r c = 0.386Mpc, and (3 = 0.705 (Mohr, Mathiesen & Evrard 1999) and 
kT = 8.5keV (Fusco-Femiano et al. 1999). We assume that n e (r) = 0.0 at r > 2Mpc. This outer boundary is roughly 
equal to that of the observed thermal X-ray emission. As a result, the total number of thermal electrons is 2.22 x 10 43 
and the total thermal energy of electrons is 4.52 x 10 63 erg. We assume B = lfiG. Evolution of non-thermal electrons 
hardly depends on the magnetic field strength unless magnetic field energy density becomes comparable to or more 
than CMB one. 

For the non-thermal electron populations, we assume that constant injection for a fiducial period t; n j and that 
momentum dependence of the injection rate is power-law. We assume that this single power-law form spectrum is 
valid between p ei mi n <p e < p e ,max, where p e ,m\n is taken to be the normalized momentum corresponding to E e = 3kT 
and Pe.max is assumed to be 10 7 , which is comparable to the critical energy of typical rich clusters up to which shock 
accelerated electron momentum spectra are power- low (Loeb & Waxman 2000; Totani & Kitayama 2000). Thus, 
Q(p e ,r,t) is, 

«*.,-.*>- {*<'*' S^f)^ CH> 

where Qq has a dimension of number density per time. Some candidates arc proposed for the origin of a non-thermal 
population. As for a electron component, however, shocks in ICM caused by cluster merger is believed to be the 
most promising because they have enough and more energy than other candidates such as AGNs or normal galaxies. 
Furthermore, it is difficult to explain the typical spatial size of radio halos and relics for point sources like AGNs and 
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galaxies because diffusion length within ICS cooling time is too short. In the present model, thus, we set ti n j = lGyr, 
which is roughly equal to the duration of shocks in cluster mergers (Takizawa 1999, 2000). 

Spatial distribution of electron injection is not clear because it depends on detailed parameters of cluster mergers. 
Thus we set fairly simple assumption. We assume it is proportional to local electron thermal energy. Total local 
injected energy of the non-thermal population is assumed to be 5% of the total local energy of thermal electrons. 
Therefor, normalization of source term Qo(r) is, 

(Mr) - o.o 5 x (Sfil) ( a^) (jTjflvSFT- 1)*.) "'■ ("J 

According to the standard theory of first-order Fermi acceleration, it is expected that exponents of momentum spectra 
fjj are 2 in non-relativistic strong shocks (see Jones & Ellison 1991). However, fairly low Mach number shocks are 
expected in cluster mergers because the ICM of the pre-merger clusters is already hot (Takizawa 1999, 2000; Miniati 
et al. 2000). In fact, some X-ray observations support this idea (e.g., Markevitch, Sarazin & Vikhlinin 1999; Kikuchi 
et al. 2000). Thus we consider the cases of [i = 2.0,2.5,3.0,3.5, and 4.0 

3. RESULTS 

3.1. Evolution of Electron Momentum Spectra 

Figure ^| shows evolution of a total electron momentum spectrum in the p, = 2.0 model. In all models electron 
momentum spectra evolve in a similar way as follows. While fresh non-thermal electrons are provided (0 < t < lGyr), 
the spectra change from the original single power-law form (oc p"^) into a steady state solution where the injection 
balances with the cooling. In the higher energy range the spectra change from N(p e ) ocp~ M into p e because ICS 
cooling is dominant. In the lower energy range, on the other hand, the spectra change into N(p e ) oc pJ AI+3 / (pf, + 1) 
because Coulomb loss is dominant. After the injection stops (t > lGyr), inverse Compton cut-off appears and shifts 
towards lower energy as time proceeds. In the lower energy range, the spectra change into a steady state solution 
without a source term, N(p e ) oc p\j (p\ + 1). 

3.2. Evolution of Hard X-ray radiation 

Figure || shows evolution of hard X-ray (20-80 keV) band total luminosity. Solid lines indicate NTB components 
and dashed lines show ICS components. The flatter source electron spectra result in more ICS and less NTB emission. 
Both components are comparable to each other in hard X-ray when /i ~ 3. Sarazin & Kcmpner (2000) shows that 
both components are comparable to each other with [i ~ 2.7 when electron momentum spectra are single power- law. 
However, note that we consider different situations In our models, electron sources are power-law (0 < t < lGyr) or not 
active (t > lGyr) and evolution of electron spectra are fully calculated. In Sarazin & Kempner (2000), on the other 
hand, power-law electron spectra are assumed for "Power-law" models, which correspond to only a very early epoch 
(t ~ 0) in our models where cooling does not play a significant role. They also consider "Cooling Electrons" model, 
where spectra are assumed to be in a steady state without any sources. This corresponds to a phase where source 
is not active (t > lGyr) in our models. Our results arc different than theirs because electron momentum spectra are 
much modified from a single power-law form in the source term because of the various cooling processes. In particular, 
spectral flattening in lower energy owing to Coulomb cooling is important. 

Although total energy injected into a non-thermal electron population is the same value in all models, the flatter the 
source electron spectra are, the more luminous the total hard X-ray emission is. This is because the electrons related 
to ICS components lose their kinetic energy mainly through ICS loss while the electrons related to NTB components 
lose their kinetic energy mainly through not bremsstrahlung loss but Coulomb interaction with the thermal ICM. 
Thus, ICS dominant models are much more efficient for energy transformation from non-thermal electrons into hard 
X-ray radiation. This point is very important when the origin of non-thermal electrons is discussed. 

After injection stops it > lGyr) , both components of hard X-ray emission decrease in all models. However, timescales 
of decline are rather different between NTB and ICS components. In all models ICS components decrease in the same 
way. The decline timescale of ICS components is determined by ICS cooling time of the relevant electrons, which is 
commonly ~ a Gyr for all models because energy of the related electrons is ~ GeV (see Fig. 0b). However, this is 
not the case for NTB components. The flatter the source spectra are, the more slowly NTB components decrease. 
The reason of this behavior is as follows. When the source spectra are flatter, equilibrium electron spectra are also 
flatter. Thus, not only weakly relativistic electrons but also higher energy (~ 100 MeV) electrons contribute to the 
hard X-ray through NTB significantly. While cooling time of weakly relativistic (~ 10 — 100 keV) electrons is very 
short (~ 10 6 ~ 7 yr, see Fig. [l]b), that of ~ 100 MeV electrons is fairly long (~ a few Gyr). As a result, the contribution 
from higher energy electrons survives for several Gyrs in NTB hard X-ray. This is more prominent when the source 
spectra are flatter. 

Observed hard X-ray spectra are often parameterized by a single power-law model. Thus, we fit our resultant total 
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Fig. 2. Evolution of a total electron momentum spectrum in the /i = 2.0 model. The different lines correspond to different epochs: 
t = 1.0 Gyr, solid line; t = 1.5 Gyr, dashed line; t = 2.0 Gyr, dot-dashed line. While fresh non-thermal electrons are provided 
(0 < t < lGyr), the spectra change from the original single power-law form into a steady state solution. In the higher energy range 
the spectra change from N(p e ) ocpj 2 into pj 3 because of ICS cooling. In the lower energy range, on the other hand, the spectra 
change into N(p e ) oc p e / (p 2 + 1) because of Coulomb loss. After the injection stops (i > lGyr), inverse Compton cut-off appears and 
shifts towards lower energy as time proceeds. In the lower energy range, the spectra change into a steady state solution without a 
source term, N(p e ) oc Pe/(Pe + 1). 

hard X-ray spectra into a single power-law model (L e oc e Q ) to get a spectral index a. Figure |3] shows evolution of the 
spectral indices of total hard X-ray (20-80 keV) spectra. The behavior depends on whether ICS or NTB components 
are dominant. When ICS components are dominant (/x < 3.0), spectral indices a changes from — (/.t — l)/2 into — /.t/2 
while the source is active (0 < t < lGyr). After the injection stops, the indices become much smaller values because 
the inverse Compton cut-off enters a hard X-ray range. When NTB components are dominant (/i > 3.0), on the other 
hand, the indices evolve in a rather different way. They become steady state values more quickly. In this stage, 
the indices are ~ —0.8 within our models. After the injection stops (t > lGyr), they become ~ —0.15, which is an 
approximation to a spectrum that varies as the logarithm of the photon energy (see Sarazin & Kempner 2000). In 
this way, NTB dominant models result in rather flat hard X-ray spectra. 

NTB flux is proportional to a volume integral of the square of the ICM density because we assume that injection 
of non-thermal electrons is proportional to the local ICM density. On the other hand, ICS flux is proportional to a 
volume integral of the ICM density itself. As a result, this ratio is proportional to density weighted mean density of 
ICM <n e > defined as follows, 

f n e (r) 2 r 2 dr , . 

<n e >^ \ e )\ 2 , ■ (17) 
J n e {r)r z ar 

X-ray observations shows that a core radius r c is correlated with central density n e fi. Fujita & Takahara (1999a, b) 
shows 
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though a dispersion of a factor of three or four is seen in the data. Using this correlation, we plot < n e > as a function 
of r c in figure (||). Solid and dashed lines indicate cases of f3 — 0.6 and 0.7, respectively. A square point shows the 
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Fig. 3. Evolution of total hard X-ray (20-80 koV) luminosity of each model. Solid and dashed lines show bremsstrahlung and inverse 
Compton components, respectively. The natter source electron spectra result in more ICS and less NTB emission. Both components 
are comparable to each other in hard X-ray when fi ~ 3. 

Coma cluster. This figure shows < n e > can become about twice of that of Coma in relatively high density clusters 
which follow the correlation. If we take a dispersion in the correlation into account, < n e > can be nearly ten times 
larger than that of Coma. Therefore, the flux ratio of NTB to ICS can be about ten times higher than the values in 
figure 3 for relatively high density clusters. 

4. COMPARISON WITH COMA CLUSTER 

In this section, we compare our models with the Coma cluster. We compare a spectral shape of hard X-ray and 
a luminosity ratio of hard X-ray to EUV of our models with Coma. Then we seek an energy injection rate which 
can reproduce Coma for each model with different values of \x and discuss their implications. In Coma, a power-law 
model gives flux of 2.2 x 10 -11 ergs cm -2 s _1 in the 20-80 keV band (Fusco-Femiano et al. 1999). This corresponds to 
luminosity of 5.1 x 10 43 erg s _1 in the same band. The spectral index of hard X-ray from Coma is poorly determined. 
0.3 > a > —1.5 is inferred at the 90% confidence level. In our parameter range of \x (2.0 <fi< 4.0), thus, all models can 
reproduce this fact while fresh non-thermal electrons are provided (see a region of < t < lGyr in Figure [I|). However, 
it is difficult for ICS dominant models to reproduce the above-mentioned range of a after the injection stops (see a 
region of t > lGyr in Figure |]). 

Excess EUV radiation is also detected in Coma. The luminosity ratio of hard X-ray to EUV can be used to constrain 
our models. In Coma, the excess EUV radiation is significantly detected only within r < 0.64Mpc and a power-law 
model gives EUV luminosity of 1.5 x 10 42 erg s" 1 in the 65-245 eV band (Bowyer, Breghofer & Korpela 1999). This 
gives a luminosity ratio I/hxr/^euv = 34. We calculate EUV luminosity due to ICS within r < 0.64 Mpc for each 
model and compare the ratios with this value. Figure ^ shows evolution of £hxr/£euv for each model (solid lines). 
An observed value of Coma derived from Fusco-Femiano et al. (1999) and Bowyer et al. (1999) is shown by dashed 
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time [Gyr] time [Gyr] 

Fig. 4. Evolution of spectral indices in hard X-ray (20-80 keV) of each model. The behavior depends on whether ICS or NTB 
components are dominant. When ICS components are dominant (/i < 3.0), spectral indices a changes from —{ji— l)/2 into —fJ./2 
while the source is active (0 < t < lGyr). After the injection stops, the indices become much smaller values because the inverse 
Compton cut-off enters a hard X-ray range. When NTB components are dominant (/i > 3.0), the indices are ~ —0.8 while sources 
are active. After the injection stops, they become ~ —0.15. 

lines. Note that £hxr/ Lew of our models must be more than the observed value because ICS component from an 
old non-thermal electron population and/or thermal emission from warm gas possibly contribute to the observed EUV 
emission, which are not considered in our models. From figure |^ it is clear that this constraint is not satisfied when 
/i = 2.5 and 3.0. When fi = 3.5, this can be satisfied only in a very early epoch (t < 0.1 Gyr). In other words, these 
models emit too much EUV compared with hard X-ray. On the other hand, models with /i = 2.0 and 4.0 satisfy the 
above-mentioned constraint in an active phase (0 < t < lGyr). The reason is as follows. When ICS is dominant in 
hard X-ray (/i < 3), the ratio of hard X-ray to EUV decreases as source spectra are steeper. This is because electrons 
relevant to hard X-ray have higher energy than those relevant to EUV. On the other hand, when NTB is dominant 
(/i > 3), that ratio increases as source spectra are steeper because electrons relevant to hard X-ray have lower energy 
than those relevant to EUV. As a result, fi = 2.0 and 4.0 models are favorable while [x = 2.5,3.0, and 3.5 models are 
not. 

Although both ICS and NTB dominant models can be consistent with the observed hard X-ray spectral index of 
Coma, our NTB dominant models do not produce observed radio spectrum. For example, Deiss et al. (1997) shows 
the radio spectral index in the frequencies lower than 1.4 GHz is -1.16, which means fj, = 2.32 in the electron source in 
an steady state. Therefor, another non-thermal electron component with flatter spectrum is necessary even if NTB is 
dominant in hard X-ray. 

Next, we seek a numerical value of a total injected energy rate -E n th which reproduce the Coma cluster hard X- 
ray luminosity at t = 1 Gyr. This is consistent with ICM temperature distribution derived from X-ray observation 
(Watanabe et al. 1999) and some numerical calculations of cluster merger (Burns et al. 1994; Ishizaka & Mineshige 
1996) which suggest that Coma cluster is in major merger which started ^0.5—1 Gyr ago. Note that E nt h is not 
very sensitive to accurate determination of an observational epoch because hard X-ray luminosity is nearly constant 
in all our models in 0.5Gyr< t < lGyr. Then, we calculate total injected energy to the non-thermal population in 1 
Gyr, -Etot.nth, and its ratio to the total electron thermal energy, £/tot,nth/-E'th- The results for each model are shown in 
table |. 

From table 1, we find that NTB dominant models require very high rates of energy injection to explain the results 
of Coma cluster. This causes violent heating of thermal ICM unless duration of the acceleration process is very short 
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Fig. 5. Density-weighted mean density < n e > as a function of r c when we consider a correlation between r c and n e o- Solid and 
dashed lines indicate cases of /3 = 0.6 and 0.7, respectively. A square point shows a value of Coma cluster. < n e > can become about 
twice of that of Coma in relatively high density clusters which follow the correlation. If we take a dispersion in the correlation into 
account, < n e > can be nearly ten times larger than that of Coma. 

(Petrosian 2001). Thus, we conclude that it is very difficult to consider NTB components significantly contribute to 
hard X-ray emission of Coma. 

On the other hand, ICS dominant models are free from the difficulties mentioned above. The required injection 
rate is ~ i0 46 ~ 47 er g s -l. This value can be realized by cluster mergers. Thus, ICS components are more likely 
as an origin of Coma hard X-ray in point of required energy. As many authors pointed, however, the flux ratio of 
synchrotron radio to hard X-ray gives us the volume averaged magnetic field strength in this case, which must coincide 
with that determined from Faraday rotation toward individual polarized radio galaxies unless spatial inhomogeneity 
plays an significant role, (e.g., Fusco-Femiano et al. 1999; Ensslin, Lieu, & Biermann 1999; Kcmpner & Sarazin 2000). 
We should look for any other causes of this discrepancy. For example, we should consider the spatial structure of 
magnetic field, inhomogeneity of magnetic field and non-thermal electrons, or spatial anti-correlation of magnetic field 
strength and non-thermal electron density, etc. Selection bias and Galactic contamination are extensively discussed 
by Goldshmidt & Rcphaeli (1993) as the possible origin of this discrepancy. 

5. CONCLUSIONS 

We calculate evolution of a non-thermal electron population in clusters of galaxies and investigate hard X-ray 
radiation due to both NTB and ICS. We calculate cases where the electron sources are power-low in momentum 
(ocpJ M ). With the parameters of Coma, ICS components dominate NTB ones in hard X-ray radiation when /j, < 3.0 
and vice versa. This dividing line may depend on thermal ICM density distribution. NTB dominant models have 
fairly flat hard X-ray spectra. 

We compare our models to the Coma cluster. Luminosity ratio of hard X-ray to EUV is consistent with Coma when 
/i = 2.0 and 4.0. When \i = 3.5, it is consistent only in a very early epoch. However, NTB dominant models require very 
high energy injection rates and cause unusual violent heating of thermal ICM unless duration of injection is very short 
(Petrosian 2001). On the other hand, ICS dominant models require much lower injection rates, which can be realized 
by cluster mergers. Thus, we conclude that ICS is more likely as the origin of Coma hard X-ray. In this case, however, 
the flux ratio of synchrotron radio to ICS hard X-ray gives us the volume averaged magnetic field strength. This value 
must coincide with that determined from Faraday rotation toward individual polarized radio galaxies unless spatial 
inhomogeneity plays an significant role. Thus, we should look for other solutions for the apparent discrepancy between 
magnetic field strength derived from the above-mentioned two methods. It may be related to spatial configurations 
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Fig. 6. Evolution of Lhxr/£euv for each model (solid lines). Observed value of Coma derived from Fusco-Femiano et al. (1999) 
and Bowyer et al. (1999) is shown by dashed lines. Note that Lhxr/^euv of our models must be more than the observed value. It is 
clear that this constraint is not satisfied when fi = 2.5 and 3.0. When /i = 3.5, this can be satisfied only in a very early epoch (t < 0.1 
Gyr). On the other hand, models with fj, = 2.0 and 4.0 satisfy the above-mentioned constraint in an active phase (0 < t < lGyr). 

of cluster magnetic field. 

Our conclusions suggest that the momentum spectrum of the electron sources in Coma is fairly flat (/u, < 2.5). 
This is consistent with radio observations. Deiss et al. (1997) shows the radio spectral index in the frequencies 
lower than 1.4 GHz is —1.16, which corresponds to \x = 2.32 in the electron source term in a steady state (Note that 
equilibrium electron spectra are steeper than source spectra because of ICS cooling). However, this implies fairly 
high Mach number shocks. According to the standard theory of shock acceleration, /j, < 2.5 means a Mach number 
more than 3.0. This value is somewhat lager than estimated values (~ 2) of merger shocks from X-ray observations of 
other clusters (e.g., Markevitch, Sarazin, & Vikhlinin 1999; Kikuchi et al. 2000). However, estimated Mach numbers 
can be underestimated because observed temperature gradients can be underestimated owing to contamination by 
foreground and/or background ICM (Takizawa 2000, Shibata et al. 2001). Thus, it is possible that such 'hidden' high 
Mach number shocks contribute to particle acceleration in Coma. Indeed, merger shocks with a Mach number more 
than 3 do exist in numerical simulations of cosmological structure formation (Miniati et al. 2000). Another possibility 
is acceleration at external accretion shocks, which can be very high Mach number shocks because accreting gas is 
much colder than cluster virial temperature (Miniati et al. 2000). 

Another possibility is reacceleration in intracluster space due to turbulent Alfven waves and/or other physical 
process (Brunetti et al. 2001; Ohno, Takizawa, & Shibata 2001). In this case an electron spectrum can be much 
modified from an original source form. As a result, situation effectively similar to our /i = 2.0 model possibly occurs 
though it depends on detailed acceleration history of non-thermal electrons. 

In general, it is believed that non-thermal energy is less than thermal one in ICM. Only ICS dominant models with 
flat source spectra can satisfy this situation among our models. Even in this case, however, the obtained non-thermal 
to thermal energy ratio is fairly high. For example, TeV gamma-ray observation of SN1006 suggests this ratio is 
a few % (Tanimori et al. 1998), whereas our fj, = 2.0 model implies this ratio is more than 10 %. However, note 
that total energy of non-thermal electrons is sensitive to a spectral shape in lower energy. We assume that source 
spectra are power-law in momentum down to E e = 3fcT. This is valid when thermal electrons in ICM are the source of 
non-thermal electrons. However, this is not the case when old non-thermal electrons which leaked from AGNs and/or 
galaxies, and/or fossil radio plasma (Ensslin & Gopal-Krishna 2001) are the source. Once electrons are accelerated 
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k>g(-E t ot,nth) 


log(Etot,nth/-Eth) 


Model 


erg s _1 


erg 




Ai = 2.0 


46.39 


62.89 


-0.77 


2.5 


47.51 


64.01 


0.35 


3.0 


49.13 


65.63 


1.97 


3.5 


49.87 


66.37 


2.71 


4.0 


50.08 


66.58 


2.93 



up to ~ 100 McV, they have a fairly long (~ a few Gyr) life time. Thus, it is possible that lower boundary energy of 
single power-law is larger and that required energy is less. 

I thank M. Hcnrikscn and K. Nakazawa for fruitful discussion and S. Shibata for continuous encouragement. I would 
also like to thank the anonymous referee for very useful comments. 
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